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We consider the out-of-equilibrium (quasi-) particle number distributions of the Higgs and W-fields during 
electroweak tachyonic preheating. We model this process by a fast quench, and perform classical real-time lattice 
simulations in the SU(2)-Higgs model in three dimensions. We discuss how to define particle numbers and effective 
energies using two-point functions in Coulomb and unitary gauge, and consider some of the associated problems. 
After an initial exponential growth in effective particle numbers, the system stabilises, allowing us to extract 
effective masses, temperatures and chemical potentials for the particles. 



1. INTRODUCTION 

In recent scenarios of electroweak baryogenesis 
|1I2| the electroweak transition is assumed to have 
taken place at low (zero) temperature shortly af- 
ter inflation, by the effective mass-squared pa- 
rameter of the Higgs field changing sign from posi- 
tive to negative ('tachyonic'). An important issue 
in this scenario is the time needed for the system 
to reach approximate thermalisation, and the re- 
sulting effective temperature. The temperature 
must be low enough and the thermalisation suffi- 
ciently rapid that sphaleron transitions that can 
wash out the generated baryon asymmetry are 
prevented. One of the main aims of this study is 
to determine the effective temperature. 

We employ the classical approximation, which 
allows us to study fields far from equilibrium non- 
perturbatively by numerical simulation. In our 
case, this approximation should be justified |3j 
since the instability resulting from the change of 
sign in the mass-squared parameter leads to ex- 
ponentially growing, and hence large, occupation 
numbers. The full results are presented in 0]. 

2. PARTICLE DISTRIBUTIONS 

The definition of particle numbers and energies 
in an interacting field theory out of equilibrium is 
not unique. Furthermore, in a non-abelian gauge 
theory the two-point correlators used to define 
them may be gauge dependent. Here, the natural 
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choice of gauge for the Higgs doublet </> is the 
unitary gauge, where it only has one non-zero real 
component = (0, /i/-\/2)"^. For the gauge fields, 
we will study the particle distribution in both the 
unitary gauge and the Coulomb gauge diAi ~ 0. 

We use a method |5lf)| where effective particle 
numbers and energies are determined selfconsis- 
tently, in analogy with the free-field quantum cor- 
relators. For the Higgs field (with Tih = h), they 
are defined as 

= {h{k)h{-k))c{nh{k)7rh{-k))c , (1) 

^H^ li^hjkyhi^h (2) 
V {h{k)h{~k))c 

where (• • ■)c is the connected two-point function 
given by {AB)c = (AB) - {A){B), and we have 
suppressed the common time coordinate t. 

The gauge field correlators {Af{k)A''j{-k)) = 

6'^''C^^{k) can be decomposed in a transverse 
and a longitudinal part, as 

C^ik) = - ^)DHk) + '^Dtik) , (3) 

and analogously for the canonical momentum 
field E'f = FjQ. In the Coulomb gauge the gauge 
potential is purely transverse, and and oj^ can 
be defined as 

n^^^^mom. u^t^s^^-y (4) 

In the unitary gauge, the transverse and lo^ 
are found to be the same as for the Coulomb 
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Figure 1. Higgs (top) and Coulomb-gauge W 
(bottom) particle distributions for different times. 



gauge, while the longitudinal occupation num- 
bers and mode energies, assuming the form oj^'^ = 
m^ff + , are given by @] 



ni^^^Dtik)Df{k), u;^ = mls^^^. (5) 

3. RESULTS 

We have performed simulations on a volume 
= 21^m]f with g = 2/3, A = 1/9, giving 
rnj^/rriy^, — 2, with lattice spacing a = 0.35m^^. 
We initialise the system according to the "Just 
the half" scheme introduced and explained in |3] . 
The Higgs field is initialised by generating clas- 
sical realisations of an ensemble reproducing the 
quantum vacuum correlators. We assume that 
the Higgs field is in the symmetric phase ((0) — 0) 
with an effective mass parameter /i^g- = /i^ > 0, 
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Figure 2. Particle numbers of the zero-modes as 
a function of time. H: Higgs; Wt- transverse W 
in unitary gauge; Wc- transverse W in Coulomb 
gauge. 



and approximately free field fluctuations 

^'^(^^'^^^^^^ = 2j{ky = (6) 

with uj{k) = ^ [p- + k^. However, we only ini- 
tialise the unstable (|fc| < |/i|) modes. The gauge 
potential A is initialised to zero, while the S-field 
is constructed to satisfy the Gauss constraint. 

We model the transition where /i^j^ goes 
through zero as a quench, in which /z^g flips its 
sign from fj? to —fi^ instantaneously. 42 indepen- 
dent realisations of the initial conditions © have 
been generated, and the subsequent time evolu- 
tion sampled for trriH = 1, 2, . . . , 12, 20, 30, 40, 
50, 100. Nearby momenta have been averaged 
within 'bins' of size A = 0.05a"^ = O.OlTSmjy 

Figure Q] shows the particle distributions for 
the Higgs and Coulomb-gauge W fields. In 
both cases we see the occupation numbers of 
the low-momentum modes increasing exponen- 
tially up to i « 6 — 8 . At this point 
the low-momentum modes start saturating, while 
the high-momentum modes rapidly become popu- 
lated. After t « 20 mjj^ the system evolves only 
very slowly, with an approximately exponential 
distribution at late times. A similar behaviour is 
seen for the W fields in the unitary gauge. The W 
distributions in unitary and Coulomb gauge start 
out very different (the unitary-gauge distribution 



initially being close to the Higgs distribution) but 
become indistinguishable after t « 40 m^^. The 
early exponential growth and subsequent slow 
evolution is illustrated in fig. [3 which shows the 
time evolution of the zero-mode occupation num- 
bers. The exponential fall-off with momentum 
means that we do not expect problems with lat- 
tice artefacts due to cutoff effects. 

Turning to the dispersion relation, we find qual- 
itatively the same behaviour for the Higgs parti- 
cles and the transverse W-modcs in the Coulomb 
and unitary gauges. For t < 20 m'j^ there is no 
sensible dispersion relation, while for t > 30 
it approaches the form lo^ = m'^ff+ck'^, with c « 1 
for the Higgs and Coulomb-gauge W particles. In 
the unitary gauge, the particle-like behaviour is 
slower to emerge, and the slope remains smaller 
than 1 (and increasing) at the latest times. For 
the longitudinal modes, with the 'inverse' defini- 
tion (jSJ of the mode energy, the dispersion rela- 
tion approaches a straight line only very slowly, 
and we do not attempt to perform a fit to the 
data. We fit the effective energies to the form 
Lul = ck^ + mlg. For the W fields, we find that 
the intercepts in the two gauges are compatible 
and change very little from i = 40 m'j^^ onwards, 
yielding a value TOos « 0.68m^r, close to the 
zero-temperature value 0.71 mn- For the Higgs 
fields, we find an effective mass which is some- 
what smaller than the zero-temperature mass, al- 
though it appears to be increasing with time. 

Taking the effective energy from the fitted dis- 
persion relation, we then fit the particle numbers 
to a Bose-Einstein distribution, rik = [exp{{LUk — 
tx) /T) - 1]-^ We only fit to the lowest 5 non-zero 
modes, for which rife > 0.5 and thus the classi- 
cal approximation may be expected to be valid. 
The best fits at the latest time (t = 100 mJi^) are 
shown in fig. 01 We find that the effective temper- 
ature T of the Coulomb-gauge W fields stabilises 
from t fa iOmJj^ onwards at fa 0.42 m/f. In 
the unitary gauge, the temperature is lower (but 
increasing) due to the smaller slope in the dis- 
persion relation. The effective Higgs temperature 
starts out higher but drops to Th ~ 0.4 m/f at 
t = 100 mjj^. In all cases, however, a large chem- 
ical potential — slightly larger even than the ef- 
fective mass — is required to describe the distri- 
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Figure 3. Higgs and W-particle numbers at t = 
100 mjj^ as a function of the effective energy ui 
taken from the fitted dispersion relation. Also 
shown are the best fits to the Bose-Einstein dis- 
tribution. The zero-modes are outside the bound- 
aries of the plot. 



bution. This would lead to an unphysical pole in 
the distribution; however, the occupation num- 
bers of the zero-modes (which are not included in 
the fits) lie considerably below the fits, indicating 
a fiattening of the distribution relative to a BE 
form at the lowest available energies. 
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